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ABSTRACT 

We calculate the strong lensing probability as a function of the image-separation Afl in TeVeS (tensor-vector- 
scalar) cosmology, which is a relativistic version of MOND (Modified Newtonian Dynamics). The lens, often 
an elliptical galaxy, is modeled by the Hernquist profile. We assume a flat cosmology with fl^ = 1 - Qa - 0.04 
and the simplest interpolating function yu(x) - min(l, x). For comparison, we recalculated the probabilities for 
lenses by Singular Isothermal Sphere (SIS) galaxy halos in LCDM with Schechter-fit velocity function. The 
amplification bias is calculated based on the magnification of the second bright image rather than the total of 
the two brighter images. Our calculations show that the Hernquist model predicts insufficient but acceptable 
probabilities in flat TeVeS cosmology compared with the results of the well defined combined sample of Cosmic 
Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS); at the same 
time, it predicts higher probabilities than SIS model in LCDM at small image separations. 
Subject headings: cosmology: theory — dark matter — galaxies: mass function — gravitational lensing — 
methods: numerical 



L INTRODUCTION 

Since Bekenstein proposed the relativistic, modified New- 
tonian dynamics (MOND) theory, named tensor-vector-scalar 
(TeVeS; Bekenstein 2004), it has become possible to inves- 
tigate the MOND phenomena in the cosmological sense. In 
particular, after determining the geometry and background 
evolution of the Universe, and calculating the deflection of 
light due to a weak gravitational field, one can test TeVeS 
and thus M OND with gravitat ional lensi ng ( Chiu, Ko & Tia^ 
120061 IZhao et al 2006: .Angus. Famaev & Zhao 2006 ). Be- 
fore TeVeS, strong gravitational lensing in the MOND regime 
could only be manipulated by extrapolating non-relativistic 
dynamics (Oin, Wu, & Zou 1995; Mortlock & Turner 2001*), 
in which th e deflection angle is only half the value in TeVeS 
dZhao & Oin.2006,). 

Needless to say, comparing the predicted results of grav- 
itational lensing with observations is of key importance in 
testing TeVeS. Zhao et al. (2006) first examined the consis- 
tency of the strong lensing predictions in the TeVeS regime 
for galaxy lenses in the CfA-Arizona Space Telescope Lens 
Survey (CASTLES). In this Letter, we investigate the statis- 
tics of strong lensing in the TeVeS regime, and compare the 
predicted lensing probabilities to the well defined sample of 
CLASS/JVAS survey. We adopt the mass function of the 
stellar component of galaxies ( P anter. Heavens. & .Timene^ 
11)04). As a first approximation, we do not consider 
galaxy cluster lenses; the lenses in the well defined sam- 
ple in CLASS/JVAS are believed to be produced by galax- 
ies rather than galaxy cllusters, although a cluster lens, 
SDSSJ10 04, was discovere d in Sloan Digital Sky Survey 
(SDSS) dlnadaetalJ l2()03t |Oguri et alj l200l . We con- 
sider the simplest MOND interpolating function jiix) and 
use the Hernquist profile (Hernauist 1991) to model the 
galaxy lenses. It is now established that, in standard cos- 
mology (LCDM), when galaxies are modeled by a Singu- 
lar Isothermal Sphere (SIS) and galaxy clusters are mod- 
eled by a Navarro-Frenk- White (NEW) profile, the pre- 
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dieted strong lensing probabilit ies can match the results of 
CLASS/JVAS quite well (e.g., "Chae" »2003^, "Chen" '2003 alH 
2004a b; Li & Ostriker 2002; Mitchell et al. 2005; Oguri et ^ 
2002; Oguri. Suto & Turner 2003; Oguri & Keeton 2004 
Pengetal. 2006; Sarbu. Rusin. & Ma 2001; Wang 20o| 
Zhang 2004; Zhang et al. 2005). For comparison, we recal- 
culate the lensing probabilities predicted by the SIS mod- 
eled galaxy lenses in LCDM cosmology with the velocity 
function. Note that, in LCDM, baryon infall effect (e.g., 
iKochanek & Whilte 2001; Keeton 2001) h as been well de- 
scribed by SIS m odel for galaxies ( Rusin & Kochanek''2003 
iKoopmans et al.l2006.) . at least statistically; furthermore, the 
effects of substructures (jOguri 2006) are also considered since 
we use the velocity function to account for the number den- 
sity of lensing galaxies. Throughout this Letter, we assume 
the source QSOs have a redshift of Zs - 1.27. 

2. TEVES COSMOLOGY AND DEFLECTION ANGLE 

As in Bekenstein (2004) and Zhao et al. (2006) we adopt 
the Friedmann-Robertson- Walker (FRW) metric in TeVeS, 
i.e., dT^ = -c^dt^ + a(tf[dx^ + f(x)(d0^ + sin^ 0di^^)] all in 
physical coordinates, where c is the speed of light and f(x) - 
X for a flat universe. The proper distance from the observer 
to an object at redshift z is DP{z) = c -i- z)LKj)Y^dz, 
where the Hubble parameter at redshift z is li(z) = d/a * 

Hq s/Hbil + z)^ + Ha, where Q.h and Oa are the constant den- 
sity parameter for baryon and dark energy, respectively, and 
we set the co ntribution from the scalar field to be zero by ap- 
proximation ^BekensteiDll2004^ ■ The angular diameter dis- 
tance from an object at redshift zi to an object at redshift Z2, 
is D{zi,Z2) = [c/(l + zj)] p Hodz/H{z) for flat universe. We 
assume cosmologies with the baryon density = 0.04 and a 
Hubble parameter /; = 0.73. 

In TeVeS, the lensing equation has the same form as in gen- 
eral relati vity (GR), and for a spherically symmetric density 
profile (iZhao et al.l2006h 

/3^0--—a, a(b)^ — — dl, (1) 

Ds Jo c^r dr 
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Fig. 1 . — Source-image relation in a flat TeVeS cosmology with z/ = 0.05 
and Zs = 1.27 for a Hemquist lens galaxy with a mass M = l.lMj, (panel a) 
and M = O.IM^ (panel b). Note that eq. (3) predicts a smaller r;,/ro hence 
a more efl'ective lens for the p = 0.35 model (soUd Une) than the p = 
model (dashed line) when M > 0.287M* (as in panel a) and vice versa when 
M < 0.287M* (as in panel b). 



where /3, 6 - bjDi and a{0) are the source position angle, 
image position angle and deflection angle, respectively; b is 
the impact parameter; D^, and Du are the angular di- 
ameter distances from the observer to the lens, to the source 
and from the lens to the source, respectively; g{r) - dO(r)/dr 
is the actual gravitational acceleration, <l)(r) is the spherical 
gravitational potential of the lensing galaxy and I is the light 
path. It is well known that the stellar component of an el- 
liptical galaxy can be well modeled by a Hernquist profile 

p(r) — 2„ru-+r, )3 ' ^''■^ '■^^ mass interior to r as M(r) - jf-^^, 
where M - 47Tr^p(r)dr is the total mass and r/, is the 
scale length. The corresponding Newtonian acceleration is 
gjv(r) = GM(r)/r2_= GM/(r + r,,f. According to MONO 
( iMilgr om 1983; ISanders & McGau gh 2002; Sanders 2006), 
g{r)fi(g{r) / ao) = gN{r)- We choose the simplest interpolating 
function fi{x) with = x for x < I and fi{x) = 1 for x > 1 . 
Thus, the deflection angle is 
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for b < ro, 

for b > ro, 
(2) 

where ro and vo are defined by GM/{ro + r/,)^ = VQ/(ro + r/,) - 
ao = 1.2 X lO^'^cmi"^, so that ro is a transition radius from 
the Newtonian to the Mondian regime, vo is the flat part of 
the circular velocity (i.e., the circular velocity in the Mon- 
dian regime). The above deflection angle has an analytical but 
cumbersome expression ( Zhao et al..2006) . so we calculate it 
numerically. 

We also need a relationship between the scale length r/, and 
the mass M, which could be determined by observational data. 
Firstly, the scale length is related to the effective (or half-light) 
radius Rg of a luminous galaxy by r/, = 7?,./ 1.8 (Hernauisn 
119911) . It has long been recognized that there exists a cor- 



relation between Rg and the mean surface brightness (/(.) in- 
terior to 7?^ (Djorgovski& Davis 1987): R^ oc (/)j0.83±o.08 
Since the luminosity interior to R^ (half-light) is - L/2 - 
K{I)eR^g, one immediately finds R^ oc L}-^^. Secondly, we need 
to know the mass-to-light ratio T = MjL oc W for elliptical 
galax ies. The observed data gives p - 0.35 (Ivan der Marell 
Il99ll): a ccordin g to MOND, however, we should find p ^ 
jSandersll2()()6l) . In any case we have L oc M^'^^^pK There- 
fore, the scale length should be related to the stellar mass 
of a galaxy by r/, oc m'-^^^^'^^'. In our actual calculations, 
we need to know ro/r/,. Since (ro + r/,) oc M'^^, we have 
'sin ^ AM-P', where p' = -0.5 + 1 .26/(1 -H p), and the coef- 
ficient A should be further determined by observational data. 
Without a well defined sample at our disposal, we use the 
galaxy lenses which have an obse rved effective radius R^ (and 
thus r ;,) in the CASTLES survey iMunoz. K ochane k. & Falcol 
fT999l) . which are listed in table 2 of lZhao et alJ J2006h . The 
fitted formulae are 



where 
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for p - 0.35, 
for = 
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7.64 X 10"Vi ^Mo is the characteristic mass of 



galaxies ( P anter. Heavens. & .Jim enez 2004). 

Figure n shows us the cases when a lens is located at red- 
shift z - 0.05 but with different values of ro/r/, and mass 
M. Here we allow /3 and 6 to take negative values due 
to symmetry. Generally, three images are produced when 
/3 < Per, where jScr is the critical source position determined by 
dpi do - and 9 <Q. For afl plausible range of /? = - 0.35, 
Figure^shows us that a smaller scale length results in a larger 
value of Per, as expected. 

3. LENSING PROBABILITY 

Usually, lensing cross section defined in the lens plane 
with image separations larger than AO is cr(> A0) - 
n:Dlpli.@[A6(M)- A9], where ©(x) is the Heaviside step func- 
tion. This is true only when A6(M) is approximately con- 
stant within Per, and the effect of the flux density ratio q,- be- 
tween the outer two brighter and fainter images can be ig- 
nored. From Figure[l]we see that this is not true, in particular 
for low mass galaxies. As usual, we consider the outer two 
images (the central image is very faint). For a given p, the left 
one (with 6* < 0) is closer to the center and is fainter, the right 
one (with 6* > 0) is further away from center and is brighter. 
On the other hand, when p increases, the fainter image ap- 
proaches to the center and becomes fainter, and opposite for 
the brighter image, so larger p corresponds to larger flux den- 
sity ratio, as is well known. We thus introduce a source posi- 
tion quantity /?^, determined by 
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0<fl<fl„. 



where 0q - 0(0) < 0, the absolute value of which is the Ein- 
stein radius, and 9cr is determined by dp/dO - for 6* < 0. 
Equation means that when yS^y,. < P < Per, the flux den- 
sity ratio would be larger than qr, which is the upper limit 
of a well defined sample. For example, in the CLASS/JVAS 
sample, qr < 10. The flux density ratio effect is strongest 
for intermediate redshift and low mass lensing galaxies, e.g., 
for z ~ 0.5 and /i(M), PgJPer ~ 0.35 at M = O.IM* and 
Pgr I Per ~ 0.15 at M = 0.01 M*. This effect can be ignored 
when the redshift of lensing galaxies z ~ or z ~ z,. On 
the other hand, we adopt the suggestion that the amplification 
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bias should be calculated based on the magnification of the 
second bright image of the t hree images rather tha n the total 
of the two brighter images ("Lopes & Millei<'2004"). For the 
source QSOs having a power-law flux distribution with slope 
f (- 2.1 in the CLASS/JVAS survey), the amplification bias 
is B(fi) = jiy-' (Oguri, 2002), where /i^S) = ||||fi„<fl<fl„. 

We thus write the lensing cross section with image- 
separation larger than A6 and flux density ratio less than 
Or and combined with the amplification bias B(J5) as 
liSchneider. Khlers & Falcnl1992tlChenl2004ah 

cr(> AO, < qr) - 

' ]f " Pii^- ' {J3)dl3, for AO < AOq, 

(if" - X^'j/^A^"' W'^A for A0O <Ae< A6,P) 



2nDi 



10, 



for AO > A0„ 



where p^e is the source position at which a lens produces 
the image separation A^, Aflo = A6'(0) is the separation of the 
two images which are just on the Einstein ring, and AOq^ - 
A6(J3q^.) is the upper-limit of the separation above which the 
flux ratio of the two images will be greater than qr- 

Now we can calculate the lensing probability with image 
separation larger than AO and flux density ratio less than qr, 
in TeVeS cosmology, for the source QSOs at mean redshift 
= 1 .27 lensed by foreground elliptical stellar galaxies by 
(e.g.. lWu.l99d) 



P(> AO, < qr) = 
dDP{z) 



— T--^dz n(M,z)(l + zfo-(> A0,< qr)dM,(6) 

Jo dz Jo 

where M,„ax is the upper limit of the mass for the lens- 
ing galaxies, and n(M,z) is the comoving number den- 
sity of galaxies for which we use the well fitted mass 
function of the stellar component of galaxies in SDSS 
given bv iPanter . Heavens. & Jimenez (2004): n(M)dM - 

«,(^)%xp(-^)^, where = (7.8 ± 0.1) x 

10-^/i^Mpc-\ a = - L159 ± 0.008 and = (7.64 ± 0.09) x 
10'"/! "^Mq. The exact value of Mmax is unimportant but we 
adopt Mn,ax = lOM* so that we do not consider the contribu- 
tion from galaxy clusters; unlike the galaxies these are domi- 
nated by gas of mass > lO'^M© with a y6-profile. 

The numerical results of equation (|6} are shown in Figure|2] 
The solid line represents the probabilities when ro/r;, = 6.3 x 
(M/M^,)"*"'' - 1 for M/L ^constant supported by MOND, 
and the dotted line represents ro/r/, - 9.5 x (M/M*)"""^^ - 1 
for M/L oc lP-^^ from observations. For comparison, the sur- 
vey results of CLASS/JVAS (Mvers et al.'2003VB rowne et all 
12003: i'aiiiaik et al. 1992: King et al. 1999) and the predicted 
probability for galaxy lensing by SIS profiles in LCDM 
are also shown. The observational probability Pobs(> A0) 
(IChen 2003b, 2004a, 2005) is plotted as a thick histogram 
in Figure |2 We recalculate the lensing probability with 
image separation larger than AO and flux density ratio less 
than qr, in flat LCDM cosmology (Q„, = 0.3 and Qa - 
0.7), for the source QSOs at mean redshift z, = 1.27 
lensed by foreground SIS modeled ga laxy halos ( Ch ae et al.l 
l2(M lMall2?)03t iMitchell et alJl2005l) : Ps,si> A0,< qr) = 
£' dz^^ l2dvn(v,z)o-sis(v,z)B, where h(v,z)dv = n*(l + 

z)^(-^)'' ^W[-(-^)^]$f- is the number density of galaxy ha- 
los at redshift z with velocity dispersion between v and v + dv 
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Fig. 2. — Predicted lens probability in TeVeS with an image separation angle 
> A6 and the flux ratio < qy = 10 by Hemquist galaxies. For comparison, 
the survey results of CLASS/JVAS (thick histogram) and the predicted prob- 
ability for lensing by SIS halos in LCDM with difl'erent amplification bias 
(B = 3.976 for dashed line and 1.09 for dash-dotted line) are also shown. 



dMitchell etal. 1 1200,51) . o-sis(v,z) = 167r^(p^( 
lensing cross section, v^g - 4.4 x lO^'^i-^) 
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is the min- 



imum velocity for lenses to produce image separation > AO 
and B is the amplification bias. We adopt (n*, vv, a',j6) - 
(0.0064/i^ lVlpc'\ 198 k ms'', -1.0,4.0) for early-type galax 



ies from IChae et all ( l2002h . According to equation ^ 
and equation (|5jl, for SIS model and CLASS sample, B - 

2 i^'" Pju(J3)^-^d/3, w ith /3q,. = 9/11 and y = 2.1 (see also 
Mitch ell et alJ2005h . Therefore, B = 1.091 based on the mag- 
nification of the fainter image (dash-dotted line in Figure |3 
and B = 3.976 based on the total magnification of two images 
(dashed line). 

4. DISCUSSION 

It has been held that, in the MOND regime, the effect of 
lensing is inefficient, in particular, that strong lensing never 
occurs (e.g., Scarpa 2006). Our calculations shown in Figure 
12 indicate, however, that this is not true. Although the Hern- 
quist model predicts insufficient lensing probabilities in a flat 
TeVeS cosmology compared with the result of CLASS/JVAS, 
the result is acceptable considering, at least, that the lensing 
galaxy can be modeled by steeper slopes and more efficient 
MOND yu-functions. 

Our results argue that TeVeS (and thus MOND) generates 
lenses with higher efficiency than CDM if the latter is mod- 
elled by SIS profile and in both cases the amplification bias 
is calculated based on the magnification of the second bright 
image (for SIS, the fainter image). Usually, B is calculated 
based on the total magnification of the two images cosidered. 
Because we introduced a cutoff /3q^. due to the flux ratio qr, the 
total magnification is 2 ~ + 1 times larger than that of the 
second bright image (depending on /3 when /3 < fiq^), which 
results in the corresponding value of B about 4 times larger 
(both for SIS in LCDM and Hernquist in TeVeS). As shown 
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in Figure 12] the lensing probabilities for SIS halos in LCDM 
with B = 3.976 (total) matches the results of CLASS/JVAS 
quite well (dashed line). Similarly, if we apply the total mag- 
nification to B for Mondian Hernquist model, the final lens- 
ing probabilities would be overpredicted compared with the 
results of CLASS/JVAS. Note that an SIS profile is more con- 
centrated in mass than a Hernquist profile, so if both pro- 
files are applied in the same regime (LCDM or TeVeS), the 
SIS profile would be more effective at lensing than Hern- 
quist. Therefore, the fact that the probabilities for SIS model 
in LCDM (B - 1.09, dash-dotted line) are lower than the 
Hernquist model in TeVeS shown in Figure |2l implies that 
MOND demonstrates a higher lensing efficiency than CDM. 
This phenomena is, in fact, not difficult to understand. It is 
well known that MOND, as an alternative to dark matter for 
solving the "missing mass" problem, takes eff'ect in the re- 
gion surrounding the luminous matter with r > ro, where a 
CDM halo is assumed to have non-zero density and its ac- 
celer ation dominates over lumin ous matter in LCDM cosmol- 
ogy ("Kaplinghat & T urneril2002l) . The deflection angle a(b) 
with impact parameter b > ro can be calculated using New- 
tonian CDM gravitation or Mondian luminous matter gravita- 
tion. We know that the acceleration g{r) in the equation for 
deflection angle for an SIS modeled CDM halo is g{r) oc r 
independent of b. So, the image separation is independent of 
the source position angle /? (when /? < p^r)^ as is well known in 
the SIS model. However, for a lensing galaxy (with no dark 
matter) modeled by a Hernquist profile, we have g(r) oc r"' 
only when r > ro (Mondian regime). So the higher probabil- 



ities indicate a higher lensing efficiency between MOND and 
CDM. 

As a first attempt at investigating strong lensing statistics 
in the TeVeS scenario, we have used the simplest interpolat- 
ing function /i(x). The deflection angle is, of course, sensitive 
to ij.{x) (Zhao & Famaey 2006). The simplest fi(x) adopted 
in this Letter corresponds to the lowest physical (or "true") 
acceleration g{r). Any offier forms of /i(x) all give stronger 
physi cal accelerations than the simplest one (Zhao&Tia^ 
"2006^. Furthermore, strong lensing is very sensitive to the 
concentration and the slope near the center of the density pro- 
file of lensing galaxies. The Hernquist model and an NFW 
model has p{r) ~ r"' near the center, both are inefficient in 
lensing. If elliptical galaxies were modeled as pressure- sup- 
ported Jaff'e model, e.g., with p(r) ~ ^^(r+ay- where a is a core 
scale length, then the lensing probability would be increased. 
Another important point is that we have assumed a flat uni- 
verse with - 0.96. However, by fitting to high-z SN la 
luminosity modulus, Zhao et al. (2006) showed that an open 
universe is more likely in TeVeS (with - 0.5). In sum- 
mary, it is promising to constrain TeVeS vs CDM through 
lensing statistics. 
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